Draft version March 2, 2013 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



THE LOW-MASS STELLAR POPULATION IN L1641: EVIDENCE FOR ENVIRONMENTAL DEPENDENCE 

OF THE STELLAR INITIAL MASS FUNCTION 

Wen-Hsin Hsu\ Lee Hartmann\ Lori Allen^, Jesus Hernandez^, S. T. Megeath^, Gregory Mosby'^, John J. 

Tobin'' and Catherine Espaillat'' 

Draft version March 2, 2013 



(N 

O 
(N 

< 



O 



> 

O 

m 

o 



ABSTRACT 

We present results from an optical photometric and spectroscopic survey of the young stellar popu- 
lation in L1641, the low-density star-forming region of the Orion A cloud south of the Orion Nebula 
Cluster (ONC). Our goal is to determine whether L1641 has a large enough low-mass population to 
make the known lack of high-mass stars a statistically-significant demonstration of environmental de- 
pendence of the upper mass stellar initial mass function (IMF) . Our spectroscopic sample consists of 
IR-excess objects selected from the Spitzer/IRAC survey and non-excess objects selected from optical 
photometry. We have spectral confirmation of 864 members, with another 98 probable members; of 
the confirmed members, 406 have infrared excesses and 458 do not. Assuming the same ratio of stars 
with and without IR excesses in the highly-extincted regions, L1641 may contain as many as ~ 1600 
stars down to ~ O.IM0, comparable within a factor of two to the the ONC. Compared to the standard 
models of the IMF, L1641 is deficient in O and early B stars to a 3-4cr significance level, assuming 
that we know of all the massive stars in L1641. 

With a forthcoming survey of the intermediate-mass stars, we will be in a better position to make 
a direct comparison with the neighboring, dense ONC, which should yield a stronger test of the 
dependence of the high-mass end of the stellar initial mass function upon environment. 



Subject headings: surveys — 
function 



stars: formation — stars: low-mass 
— stars: pre-main sequence 



stars: luminosity function, mass 



1. INTRODUCTION 

Do the most massive stars form preferentially in dense, 
massive environments? Answering this question has 
proved to be frustratingly difficult. The recent review 
by Bastian et al. (2010) concludes that there is no 
clear evidence for nonstandard stellar initial mass func- 
tions (IMFs) in specific environments, though this issue 
"clearly warrant (s) further study". One problem is that 
most of the studies of IMFs in relatively unevolved young 
regions have been conducted on young star clusters (Bas- 
tian et al. 2010; Weidner, Kroupa, & Bonncll 2010, and 
references therein), because it is much more difficult to 
study low-density, dispersed regions. Perhaps the most 
extensive survey to dat e on star forma t ion in low-density 
enviroments is that of iLuhman et al.l (|2009l ). which ar- 
gued that samples of roughly 150 and 300 members of 
Taurus (depending upon the region surveyed) showed 
an anomalous IMF, particularly in suggesting a peak 
near spectral types K7-M0 not observed in other regions. 
However, the low-mass population of Taurus is insuffi- 
ciently numerous to make a strong test of the high-mass 
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end of the IMF. 

The lower part of the Orion A molecular cloud, south 
of the Orion Nebula Cluster (ONC) region, appears to 
be one site where a significant test of the dependence 
of the upper mass distribution in a lower-density envi- 
ronment can be carried out. For simplicity we call this 
region "LI 641" , though this is a broader use of the origi- 
nal Lyn ds cloud designation (a s discussed further in § 14.21 
also see lAllen fc Davis! (|20080 for a discussion of the in- 
dividual clouds). As the lower part of the Orion A cloud 
is contiguous with the ONC region, it can be assumed 
that L1641 is at roughly the s a me distance as t he ONC 
(~ 420pc; IMenten et atl [20071: IKim et al.l [20081 ) . allow- 
ing the stellar populations of the two regions to be com- 
pared directly. Based on previous surveys, as well as the 
recent Spitzer/IRAC survey of infrared-excess stars (see 
Figure [T] for locations of the IR excess objects, data to 
be published in Megeath et al. 2012), it appears that 
L1641 contains a much larger pre-main sequence pop- 
ulation than other ne arby low-density regions such as 
Taur us- Auriga (see, e.g lLuhman et al.ll2009l : lRebull et all 

[miD. 

The star of earliest spectral type spatially associated 
with the mole cular gas south of -6.1° declination is B4 
(|Racineiri968[ ). with only an additional 12 late B stars 
proje cted on the cloud according to current surveys (jSkifj 
120 101) . The ONC, on the other hand, has two O stars 
and seven stars with spectral types earlier than B4, and 
its IMF is consistent with a Salpet er upper-mass slope 
(jMuench et al.l [20021: IDa Rio et ail [l009). The ques- 
tion to be addressed is whether this apparent deficit 
of high-mass stars in L1641 is the result of its being 
a relatively low-density environment, in contrast to the 
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Figure 1. Blue boundary lines represent the fields covered by the 
Spitzer/IRAC survey of Orion (Megeath et al. 2012) and the red 
stars show IR-oxcess objec ts identified by th is survey, which are 
overlaid on ^'^CO map from lBallv et all l|1987l ). The physical scale 
of 5 pc is also shown in the plot, assuming a distance of 420pc. 

densely clustered ONC, or whether L1641 simply does 
not have enough members in comparison with the ONC 
(with ^ 2000 members; Hillenbrand &: Hartmann 1998, 
Muench et a"Lll2002l . iDa Rio et aT]|2009l . iRobberto et~aLl 
201^ that high mass stars are probable, assuming a uni- 
versal IMF. 

In this paper, we present optical V and I band photom- 
etry and low resolution optical spectroscopy of low mass 
stars in L1641. We estimate that there are about 1600 
low-mass members in the region, assuming that the disk 
fraction is constant throughout the cloud. Furthermore, 
we show that the average age of the stars in our sam- 
ple is only slightly greater than that of the ONC region. 
Thus, our findings support the idea that the comparison 
between L1641 and the ONC ultimately will provide a 
good test of the possible environmental dependence of 
the upper mass IMF. For this paper, we consider O & 
B stars found in the literature and assume that there 
are no deeply embedded early B stars that we do not 
already know. Under this assumption, we then make 
comparisons with standard forms of the IMF which in- 
dicate that L1641 is deficient in early B and O stars, 
based on the known low-mass population. In a forth- 
coming paper, we will discuss our search for high-mass 
members and add intermediate-mass members verified by 
spectroscopy, which are currently missing from our sam- 
ple. In §2 we describe our observational program and 
the procedures used to reduce the data. In §3 we present 
an analysis of our photometric and spectroscopic results. 
Then in §4 we summarize our knowledge of the low-mass 
population, briefly relate our results to the known popu- 
lation of early-type stars, consider some uncertainties re- 
lated to deciding where to divide the region between low- 
and high-density areas, and make some preliminary tests 
using standard models of the IMF. Finally, we present 



Figure 2. OSMOS optical photometry fields (blue circles) over- 
laid on ^■^CO map l|Ballv et aLIHQSTI VThe physical scale of 5 pc is 
also shown in the plot, assuming a distance of 420pc. 

our conclusions in §5. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Optical Photometry 

We obtained V and I band optical photometry with 
the Ohio State Mufri-Object Spectro graph (OSMOS 
on the MDM 2.4m Hiltner telescope (jston et al.l 1201 
M artini et al. 2011). The detector has a 18.5' x 18.5' field 
of view, vignetted by a 20' diameter circle. We used 2x2 
binning, which gives a plate scale of 0.55" per pixel. The 
majority of L1641 from -6° to -9.2° was covered with 41 
fields (see Figure [2]) . The exposure sequence consists of 
one short exposure (5 seconds) and two long exposures 
(60 seconds) in both V and I bands. The photometry 
was obtained on the nights of Dec 7 - Dec 8 and Dec 8 - 
Dec 9, 2010. Most fields were observed with the exposure 
sequence once on each night, but a few fields were only 
observed on one night. The conditions were not com- 
pletely photometric, with some visible cirrus at the be- 
ginning of each night, and the seeing ranged from 1.1" to 
1.8" FWHM. However, we started observing L1641 three 
to four hours after twilight ended and most of the cir- 
rus had disappeared by then. We observed two Landolt 
standard fields every two hours throughout the night for 
photometric calibration. SA92, SA95, SA98 and SAlOl 
were used according to the time of the night. 

Each CCD frame was first corrected by overscan us- 
ing the IDL program proc4k written by Jason East- 
man. Note that the original program is written for the 
Ohio State 4k CCD imager alone, and has been mod- 
ified to process OSMOS data. We then performed the 
basic reduction following the standard procedure using 
IRAF. Note that we used sky flats in the flat-field cor- 
rection. We determined the astrometric solution with 
imwcs in WCSTools, using the coordinates of stars from 
the 2MASS catalog. The astrometric solutions are in 
general good within one pixel, or 0.55". We then ob- 
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Figure 3. Errors in V (left) and I (right) band photometry. The errors shown here are the greater of the IRAF calculated error and the 
median absolute deviation. 



tained aperture photometry with the IRAF phot pack- 
age, using an aperture of 8, or ~ 3-4 times the FWHM. 
We used the Landolt fields for photometric calibrations. 
The rms departures of the standard stars from the cal- 
ibration equations are ~ 0.05 mag for both V and I 
bands. The MDM 4k imager used by OSMOS has a 
known issue of crosstalk between the four CCD segments. 
When a CCD pixel is saturated, it creates spurious point- 
sources in corresponding pixels on all other three seg- 
ments. Therefore, instead of finding point sources in 
the images with daofind-like packages, we used the po- 
sitions of all 2MASS point sources for our photometry. 
Non-detections and saturated stars in either bands are 
removed. If a star is saturated in the long exposure but 
not in the short exposure, the short exposure is used 
for the photometry. Measurements that are affected by 
bleeding or crosstalk from a saturated star are manually 
removed. Finally, we combine all the photometry mea- 
surements for the same star by taking the median value. 

To estimate the photometric error, we use both the 
error estimated by the IRAF phot and the variation 
in separate observations of the same star. The IRAF 
phot package only takes into account of Poisson statistics 
which tends to underestimate the true uncertainty. To 
provide a better estimate of our errors, we take the stars 
that we have at least four long-exposure measurements 
to find the median absolute deviation (MAD), a robust 
measure of the deviation. Since most fields arc observed 
on both nights and the measurements from both nights 
are used, the deviation implicitly takes into account of 
the potential effects due to variable atmospheric extinc- 
tion and seeing from one night to another. If a star is 
saturated in the long exposures, then the short-exposure 
measurements (in the case where there are more than 
two short exposures) are used. We estimate the photo- 
metric error for each star from the greater of the phot 
calculated error and the MAD. The results are shown 
in Figure [3l The typical errors are less than 0.05 mag- 
nitude for stars brighter than V=19orI = 18. This 
error is comparable to the errors obtained in the cali- 
bration of Landolt stars and is most likely due to the 
non-photometric weather conditions. The typical error 
then increases drastically with fainter stars due to un- 
certainties in photon statistics. A small fraction of stars 



have errors much larger than the typical errors, which 
can be due to non-photometric conditions (cirrus in the 
beginning of the night) as well as intrinsic variability of 
the stars from one night to another. We think intrinsic 
variability contributes to the measured errors because 
we examined the list of stars with the largest errors and 
many of them turned out to be known variable stars. 

We were able to obtain photometry of stars with V 
and I magnitude between 12 and 22 mags. However, 
from Figure |3l we know that the typical error of a V = 21 
mag star is about 0.2 mag and the error increases rapidly 
for fainter objects. We therefore limit our photometric 
sample to stars with V magnitudes between 12.5 and 21, 
I magnitudes between 12 and 20.5 and V-I between 0.5 
and 4 mags, with the typical error of a 21 magnitude 
star being 0.2 mags. 

Figure [3] shows the V vs. V-I color-magnitude diagram 
of the sample mentioned above. The top panel shows all 
the stars in our sample, regardless of its IR-excess and 
the bottom panel shows the same plot, with the class I 
protostars shown in blue diamonds and the class II's (IR 
excess objects with disks) shown in red circles. 

We note that most of the class II stars in Figure |4] fall 
in a small region on the CMD, which we approximate 
with the top two parallel lines. The bottom of the YSO 
region is traced by a line through the two points (V-I, V) 
— (0.5, 11.5) and (4, 21.5), and the top is 2 magnitudes 
above it in V at a given V-I. We expect that the Class 
Ill's (pre-main sequence objects with no IR excess) have 
similar ages and distances and therefore fall in the same 
region of the CMD. However, not all the non-IR excess 
objects in this region are members. We therefore target 
them in our optical spectroscopy surveys to confirm their 
membership. 

2.2. Optical Spectra 

There are two parts to our optical spectroscopic sam- 
ple. The first part is the sample of infrared-excess stars 
from the Spitzer survey of Megeath et al. (2012). Fig- 
ure [1] shows fields covered by the Spitzer survey and the 
positions of the IR-excess stars. The Spitzer survey iden- 
tified 166 protostars and 723 disk objects in L1641 (de- 
fined as the field outlined in Figure [T] but stops at -6°.) 
Note that only 8% of the protostars and 48% of the disk 
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Figure 4. Top; V vs. V-I color-magnitude diagram of objects 
in L1641, showing all the photometric data. Bottom: Same figure, 
with protostars (blue diamonds) and disk objects (red circles) over- 
plotted. The YSOs with disks (red circles) lie mostly in a small 
region of the CMD, marking the region populated by the pre-main 
sequence stars. There arc 772 non-IR excess sources (grey crosses) 
in the same region of the CMD. This photometry was taken in Dec 
2010 with OSMOS on the MDM 2.4m telescope. T he ext i nction 
vector is the standard extinction taken from Cardel li et al.l l|1989l ) 
with =3.1. The parallel lines are discussed in S 14.11 

objects are found in our optical photometric sample and 
are therefore sufhciently bright for optical spectroscopy. 
The second part of our sample is non-IR excess objects 
selected based on their optical photometry. There are 
a total of 772 non-excess stars that satisfy the optical 
photometry selection and we were able to target 95% of 
these objects in our spectroscopy. The sample is, how- 
ever, biased against heavily-reddened objects. 

The spectra were obtained with multi-aperture spec- 
trographs covering extended fields of view. Approxi- 
mately 75% of our spectra were obtained with Hectospec, 
the fiber- fed mul tiobject spectrograph on the 6.5m MMT 
on Mt. Hopkins (jFabricant et al.(l2cio5j ). It has 300 fibers 
that can be placed within a 1° diameter field. Each Hec- 
tospec fiber subtends 1.5" on the sky. The observations 
were taken with 270 line mm~^ grating, providing ^ 6. 2 A 
resolution and spectral coverage of 3650 to 9200A. With 
17 observations from 2006 to 2011, a total of 2369 stars 
were observed. The total on-target exposure times are 
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Figure 5. Optical spectroscopy fields overlaid on ^'^CO map from 
Bally et al., (19873 . The Hectospec fields (top) are chosen to maxi- 
mize number of IR excess objects observed. The IMACS fields (bot- 
tom) sample are either IR-selected (blue solid circles) or optical- 
selected (magenta dashed circles). The IMACS fields are mainly 
used to compliment the Hectospec fields in crowded regions. The 
physical scale of 5 pc is also shown in the plot, assuming a distance 
of 420pc. 

45 to 50 minutes for most fields, with the exception that 
the Jan 30, 2009 field was observed for 22.5 minutes. 
The top panel of Figure [5] shows fields of the Hectospec 
observations. Targets for observing runs from 2006 to 
2009 are selected based on their IR colors (fields shown 
in blue) and the targets of the 2011 run are selected based 
on their optical photometry (fields shown in magenta, see 
ii l2.1l for selection criteria). We first maximized the num- 
ber of targets that satisfy either the IR excess or optical 
selection criteria, and then filled in the remaining fiber 
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positions with other stars in the field. 

The Hectospec data were reduced throug h the stan- 
dard Hectospec data reduction pipeline (iMink et alj 
120071 ). except for sky background subtraction. The 
pipeline assumes that the sky background does not vary 
significantly with position on the sky and uses "sky 
fibers", or fibers that point to empty portions of the 
sky, to correct for the sky background. However, in a 
star-forming region, nebulosity may result in significantly 
varying sky background and therefore bad estimates of 
nebular lines. Our sky subtraction takes into account the 
spatial variations of the sky background. Some of the 
observations were taken along with offset sky spectra, 5" 
apart from the star, in which case we subtracted the off- 
set sky spectra from the science target spectra. When 
the offset sky spectra were not taken, we used the closest 
3 to 5 sky fibers as the average sky and subtract it from 
the science spectra. We tested our sky subtraction meth- 
ods with sky fibers (so a perfect sky subtraction would 
leave a spectrum with only noises around 0) and found 
that the quality of the resulting spectra is significantly 
improved with the offset sky spectra, with most of the Ha 
emission and other sky lines accounted for whereas the 
nearest 3-5 sky fibers improves the overall sky subtrac- 
tion (especially night-sky emission lines in the near-IR) 
but does not correct for the nebular lines very well. 

Another 25% of our spectra were obtained with 
IMACS, the Inamori-Magellan Areal Camera & 
Spectrograph on the Magellan Baadc telescope 
(|Bigelow fc Dressleii [200l . We used IMACS f/2 
camera in multi-slit spectroscopy mode with the 300 
line grism at a blaze angle of 17.5°. With a 0.6" slit, 
this configuration yields a resolution of 4A and spectral 
coverage of approximately 4000 to 9000A (stars close to 
the edge of the fields may not have full spectral cover- 
age). From 17 observations in 2010 and 2011, a total 
of 715 stars were observed. The standard observation 
time for each field is 5 x 10 minutes, but we increased 
the time to 6 x 10 minutes for a few observations at 
higher airmasses. The bottom panel of Figure [5] shows 
the fields of the IMACS observations. The targets of 
the first IMACS run were selected based on their IR 
colors (shown in blue), and the targets of the second 
run are selected based on their optical photometry 
(shown in magenta, see § 12.11 for selection criteria) . The 
IMACS data were reduced with COSMOS, the Carnegie 
Observatories System for MultiObject Spectroscopy, 
following the standard cookbook. 

Under good observing conditions, we were able to 
spectral-type stars as faint as V ~ 21. However, the 
faintest magnitude we can spectral type depends on the 
seeing and the instrument used. 

3. RESULTS 

3.1. Photometry 

Figure |4] shows the V vs. V ~ I color-magnitude dia- 
gram of the stars between V — I ^ 0.5 — 4 and V < 21. 
There are a total of 4475 stars in this photometric sam- 
ple. Correlating th e optical photometry a nd the IR- 
excess classification ()Gutermuth et all I2009D , we mark 
the IR-excess objects in our photometry (bottom panel). 
13 of the stars were protostars (blue diamonds) and 347 
of them were YSOs with disks (red circles). Compared to 



the number of IR-excess objects identified in the Spitzer 
survey, only 8% of the protostars and 48% of the disk 
objects are found in our photometric sample. The op- 
tical photometry fields cover slightly less area than the 
Spitzer survey, but cover the 95% of the IR-excess ob- 
jects. It is not surprising that we do not find most of 
the protostars as they are highly embedded in their en- 
velopes. We were only able to identify half of the stars 
with disk excesses in the optical photometry due to ex- 
tinction (see 13.41) . 

3.2. Spectra and Spectral Types 

Figure[6]and[7]show some example spectra of confirmed 
members from Hectospec and IMACS, respectively, ar- 
ranged by t heir spectral types. Th e Ha equivalent width 
criteria by IWhite fc Basril (|2003D is used to determine 
whether they are classical T Tauri stars (CTTS) or weak- 
hue T Tauri stars (WTTS). Note that IMACS has very 
little transmission in the blue whereas the transmission of 
Hectospec is more uniform across the spectral range. The 
[O I] hue at 5577A is very prominent in some Hectospec 
data but is actually a sky line that was not correctly 
subtracted. In this case, the Ha emission is also likely 
affected by bad sky subtraction. Figure [5] show IMACS 
spectra of rapidly- accreting stars. They not only have 
a high Ha equivalent width, but also have other emis- 
sion lines such as [O I], He I and the infrared Ca triplet. 
The last two spectra have particularly high Ha equiv- 
alent widths and show O I emission at 7773A and the 
higher orders of the Paschen series. 

We determine the spectral types of our Hectospec 
and IMACS data with SPTCLASS, a semi-au tomatic 
spectral- typing program ([Hernandez et al.]|2004D . It uses 
empirical relations of spectral type and equivalent widths 
to classify stars. It has three schemes optimized for dif- 
ferent mass ranges (K5 or later, late F to early K and 
F5 or earlier), which use different sets of lines. The user 
has to manually choose the best scheme for each star 
based on the prominent features in the spectrum and 
the consistency of several indicators. While SPTCLASS 
is insensitive to reddening and S/N of the spectra (as 
long as one can obtain a good flux estimate), it does 
not take into account the effect of the hot continuum 
emission produced by the accretion shocks. This contin- 
uum emission makes the photospheric absorption lines 
appear weaker. SPTCLASS generally assigns an earlier 
spectral type to veiled stars than their and therefore the 
SPTCLASS outputs should be considered as the earliest 
spectral type limits. Highly veiled stars (such as the ones 
shown in Figure [8]) are not spectral-typed. Eight of our 
program stars are too veiled for spectral-typing. 

We also complement our results with data from Fang 
et al. (2009, hereafter F09), which include spectral 
type estimates and Ha and Li I equivalent widths for 
266 stars. The majority of the spectral types come 
from their VLT/VIMOS o bserva tions and som e are from 
'Galfa lk fc OlofssonI (l2008l) and lAUenl (119951) . We ob- 
served 206 of the stars in the F09 sample; Figure [HI com- 
pares our spectral types with theirs. Since the spectral 
types in F09 come from multiple sources, we use differ- 
ent symbols to indicate the source. Most of the spectral 
type are consistent within 2 subclasses and there is no 
systematic trend in the differences; the only major large 
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Figure 6. Examples of spectra from Hectospec with a zoom-in view around Ho (A 6563) and Li I (A 6707). Prom top to bottom, the 
spectra are spectral typed: K3, K6, MO, MS and M5. The [O /] line at 5577A is a prominent s ky line. The red end of the M type spectra 
are affected by series of sky lines. Their TTS type are classified using the definition given bv lWhite &: Basnl II2003I ). The Ho equivalent 
width in the K6 star is bordering the cutoff in the definition. 



inconsistency is for a rapidly-accreting G type star, which 
has an uncertain spectral type because of very high veil- 
ing. For consistency we therefore use the spectral type 
from our observations for the objects in common, and 
incorporate the data from F09 for the stars we did not 
observe. 

Figure [TU] shows the distribution of spectral types of 
confirmed members. The distribution peaks around M4; 
we are systematically incomplete in the later-type popu- 
lation because of extinction (see §3.4). Due to the mag- 
nitude criteria in our target selection, there are very few 



stars earlier than K in our sample. We will address this 
in a forthcoming paper, where we attempt to identify the 
intermediate-mass members. 

We use the following criteria to identify young stars in 
our sample: 

(1) We include all IR excess stars using Spitzer/IRAC 
and 2MASS J HK colors (the IR-exces s selection criteria 
can be found in lGutermuth et "aI1l2009l : the complete cat- 
alog will be published in Megeath et al. 2012). For all 
non-IR excess stars, we identify them as members if (2) 
Li absorption at A6707A is clearly detectable, and has an 
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Figure 7. Examples of spectra from IMACS with a zoom-in view around Ha (A 6563) and Li I (A 6707). Prom top to bottom, the spectra 
are spectral typed: K3, K6, MO, MS and M5. 



Table 1 
Spectroscopic Members 



SoTiree 


( 'oiiiiriiiecl AleiiilxTs 


IR excess 


iiti IH t!xeess 


Prt)ljaljle 


Total 


Hectospec 


517 


235 


282 


69 


586 


IMACS 


311 


142 


169 


24 


335 


FangOQ* 


36 


29 


7 


5 


41 


I'otal 


864 


406 


458 


98 


962 



" Fang et al. (2009) cataloged 266 stars in L1641. The numbers here show only 
stars we did not observe. 




Figure 8. Example IMACS spectra of rapid accretors with a zoom-in view around Ha (A 6563) and Li I (A 6707). The spectra are 
highly-veiled and therefore makes spectral typing difficult. The main emission lines have been marked. 
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Figure 9. Com parison of the spectral types obtained by 
IFang et al.l l|2009l ) and this work for the 166 objects that are in- 
cluded in both samples. The catalog of Fang et al. (2009) comes 
from three sources: 1. their VLT/VIMOS observations (crosses) 
, 2. Galfalk&Olofsson(2008) (squares) and 3. Allen (1995) (di- 
amonds). The highly discrepant star near the bottom is heavily 
veiled (see text) 



0.5 



0.0 



KO 

spectnal type 



Figure 10. Distribution of spectral types of all confirmed mem- 
bers with spectral type information from Table [2] Due to the 
magnitude criteria in our target selection, there are missing stars 
earlier than K in our sample. 

equivalent width consistent with youth (the EW(Li) has 
to be greater than the Pleiades val ue for the same spec- 
tral type from IBriceno "eFall (IT997I) '). or (3) Ha emission 
is clearly above the background nebulosity. 

Criterion (2) is reliable in selecting pre-main sequence 
M and late K stars, but for early K and G stars it can 
be more problematic, as Li depleti on timescales can b e 
quite long (see, e.g., discussion in IBriceno et aIlll997D . 
Thus, to identify G and early K members, we require 
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Figure 11. Top: Positions of spectroscopically confirmed mem- 
bers in Table[2]overplotted on ^^CO map from Bally et al. (1987). 
The symbols indicate the sources of the spectra. Red stars and 
blue circles represent objects observed with Hectospcc and IMACS, 
respectively. Green triangles represent spectral types taken from 
F09. Bottom: Positions of probable members in Table [S] 

that the equivalent widths of Li absorpti on lines to be 
great er than that observed in Pleiades (jBriceno et al.l 
|1997[ ) or the Ha equivalent widths in emission or filled 
in above the expec ted absorption equivalent widths 
(IStauffer et al.l[T997l) . 

Whether or not the observed Ha emission - criterion 
(3) - is actually from the star or simply nebular emis- 
sion can be difficult to distinguish in some regions and 
especially for faint stars. We therefore compiled two lists 
of objects, one with only confirmed members of L1641 
and one that are potential members but have a higher 
contamination rate. 
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To ensure a self-consistent sample, we extracted mem- 
bers from the F09 sample using the same criteria used in 
our own spectra; therefore not all the members listed in 
F09 are included in our results here. In total, we adopted 
the F09 data for 41 stars that were not targeted in our 
observations, five of them are considered "probable mem- 
bers" due to the lack of Li absorption data. 

Table[l]lists the number of confirmed members (includ- 
ing IR-excess members and non-excess members) and 
probable members from the three sources (Hectospec, 
IMACS, F09). Table [1 gives the positions, photome- 
try, spectral types, etc. of the 864 confirmed members, 
confirmed either by their IR excess or unambiguous Li 
absorption. Out of the 864 confirmed members, 406 of 
them have IR excess, while 458 have no excess. Eight 
of the confirmed members are highly veiled by excess 
continuum emission, thought to be prod uced by the ac - 
cretion shock on the stellar photosphere (jKoenig]||199lD . 
Spectral typing is not possible for these objects. 

Table[3]gives the positions, photometry, spectral types, 
etc. of the 98 probable members. The probable members 
are stars that show Ha emission and are in the strip of V 
vs. V-I CMD (see bottom panel of Figure that most 
of the confirmed YSOs reside in. They do not, however, 
have IR-excess or unambiguous Li absorption line (most 
likely due to low S/N of the spectra). We therefore cata- 
log them as probable members. These probable members 
are excluded in our analysis and statistics unless other- 
wise noted. 

The top panel of Figure [11] shows the positions of 
spectroscopicall y confirmed r nembe rs overplotted on the 
13 CO map from lBallv et all (|1987D . The distribution of 
the spectroscopic members in general follow the gas dis- 
tribution. The bottom panel shows the positions of prob- 
able members. Probable members that are located near 
the group of confirmed members are more likely members 
whereas the ones far away from any confirmed members 
are likely non-members with nebular Ha emission. 

Figure fT2l shows the spectroscopically confirmed mem- 
bers on the CMD color-coded by their spectral types. It 
shows that we can spectral type objects down to M5- 
M6 as long as they are not highly extincted ( in gen- 
eral. Ay < 2). There are some objects well below the 
YSO regions and most of which are IR-excess objects. 
Most likely they are edge on disk objects that are seen 
in scattered light and therefore appear much bluer than 
typical YSOs. This figure also shows that the magnitude 
limits of our optical photometry and spectroscopy are 
well-matched. 

T he solid lin e s in F igure [T^ from top to bottom, are 
the iSiess et"an (2000) isochrones for 1, 2, 3 and 4 Myr. 
Because the isochrones are almost parallel to the red- 
dening vector, we can use the non-extinction corrected 
CMD to estimate the age of the population. Most of the 
population lies within the 1 Myr and 4 Myr isochrone, 
and the distribut ion peaks arou i id 3 - 4 Myrs. Using the 
Siess isochrones, iDa Rio et al.l ()2010f ) found the age of 
the ONC to be - 2 - 3 Myr. The L1641 population, on 
average, thus appears to be slightly older than the ONC 
population, even though the age difference is compara- 
ble to the age uncertainties. This means that ultimately 
we should be able to make a direct comparison for IMFs 
using spectral types without correction. 




Figure 12. Spectroscopically confirmed members on the CMD 
color-coded by their spectral types (not corrected for extinction^ 
The s olid lines are 1, 2, 3 and 4 Myr isochrones from Sicss et id] 
1120001) (from top to bottom) and the dotted line is the line used 
to define the YSO region. The minimum mass of the plotted 
isochrones is O.IMq. 

3.3. Ha as a Membership Indicator and Accretion 
Diagnostics 

The optical spectra also provide us with diagnostics of 
accretion. Table [2] and [3] list equivalent widths of Ha, 
which can be used as an indicator of accretion as the 
Hg lumino s ity co rrelates with the accretion luminosity. 
iFang et all (|2009l ) has an extensive discussion on accre- 
tion rates and disk properties in selected areas of L1641 
using Ha, H^ and He I A 5876. 

Table |4] lists other accretion indicator emission lines we 
found in the spectra. Note that if an emission line is not 
listed, it does not necessarily mean there is no emission 
line, but rather we cannot confidently tell that there is 
an emission line due to sky subtraction issues or limited 
spectral coverage in some IMACS data. 

Figure [13] shows Ha equivalent widths of all confirmed 
members vs. spectral type. The Ha equivalent widths 
can be used as an indicator of accretion as the Ha lumi- 
nosity correlates with the accretion luminosity. In gen- 
eral, the Ha equivalent widths have an increasing trend 
toward later spectral types and the IR-excess objects 
(red open diamonds) have higher Ha equivalent widths 
than the non-excess objects (black solid triangles). The 
shaded region shows t he Ha equivalent widt h criteria for 
CTTS and WTTS bv lWhite & Basrl (pOOl . Using this 
T Tauri criteria, 299 of the confirmed are CTTS and 565 
are WTTS, which gives a CTTS/ WTTS ratio of 53%. 
In the early K to G spectral type range, there are a few 
stars with very low EW(Ha) equivalent widths. 

We now look at the stars with very low Ha equiva- 
lent widths in more detail. The top panel of Figure [14] 
shows the Ha equivalent widths vs. spectral type for ob- 
jects with Ha equivalent widths less than lA in emission. 
Note that in this plot emission is positive and that the 
y axis is in linear scale. The dashed line show the Ha 
equivalent widths mea sured in the Plei ades from chro- 
mospheric activity (St auffer et al.l [19971) . We note that 
four of the members fall on or below the dashed line. 
The bottom panel of Figure [14] shows the Li equivalent 
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Table 2 

Confirmed Members in L1641. This tabic is available only on-line as a machine-readable table. 



ID 


RA 


Dec 


V 


I 


J" 


H" 


Ks° Sp 


ectral Type 


IR. exccss^^ 


EW(IIc() 


EW(Li) 


Date observed* 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 












1 


83.21121 


-6.18593 


17.430 


14.793 


13.450 


12.740 


12.537 


M3.0±0.6'' 


N 


-4.30 


0.4 


Oct 22, 2011 


2 


83.22298 


-6.10026 


17.670 


14.743 


13.146 


12.499 


12.259 


M3.6±0.6° 


N 


-11.0 


0.4 


Oct 22, 2011 


3 


83.23454 


-6.05524 


16.566 


14.356 


13.040 


12.405 


12.235 


M1.5±0.7° 


N 


-6.00 


0.5 


Oct 22, 2011 


4 


83.24131 


-6.04519 


16.819 


14.377 


12.975 


12.324 


12.117 


M2.6±0.5° 


N 


-4.00 


0.3 


Oct 22, 2011 


5 


83.30244 


-6.05771 


17.607 


14.925 


13.508 


12.879 


12.632 


M3.0±0.6'^ 


N 


-4.20 


0.3 


Oct 22, 2011 


6 


83.30896 


-6.19689 


17.006 


14.462 


13.006 


12.329 


12.074 


M2.6±0.5'^ 


N 


-22.6 


0.3 


Oct 22, 2011 


7 


83.31703 


-6.30980 


15.290 


13.320 


12.061 


11.383 


11.151 


M0.5±0.5° 


N 


-3.60 


0.5 


Oct 22, 2011 


8 


83.33062 


-6.24046 


17.858 


14.796 


13.158 


12.535 


12.274 


M4.6±0.6'* 


N 


-7.90 


0.3 


Oct 11, 2011 


9 


83.33256 


-6.07309 


16.396 


14.081 


12.722 


12.013 


11.826 


M2.2±0.6" 


Y 


-5.90 


0.4 


Oct 22, 2011 


10 


83.34235 


-6.19836 


15.604 


13.608 


12.415 


11.724 


11.530 


M0.6±0.5'* 


N 


-3.20 


0.5 


Oct 11, 2011 


" J,H,Ks photometry is from 2MASS 






















" Criteria for IR 


excess are 


defined in 


Gutermuth et al. 


(2009) and Megeath et al. (2012) 













^ Data from Hcctospcc spectra. 
Data from IMACS spectra. 

Data from .Fang et al.i (ii2009i) . No spectral-typing error estimate. 
* Local date at the beginning of the night. 



Table 3 

Probable Members in L1641. This table is available only on-line as a machine-readable table. 



ID 


RA 
(J2000) 


Dec 
(J2000) 


V 
(mag) 


I 

(mag) 


J" 
(mag) 


(mag) 


Ks" Sp 
(mag) 


ectral Type 


IR excess^ 


EW(Ha) 

(A) 


EW(Li) 

(A) 


Date observed* 


1 


83.31244 


-6.19215 


18.897 


15.697 


14.001 


13.429 


13.142 


M4.6±0.7'* 


N 


-9.40 




Oct 11, 2011 


2 


83.33856 


-6.01960 


19.102 


15.608 


13.731 


13.171 


12.890 


M4.6±1.0° 


N 


-11.6 




Oct 22, 2011 


3 


83.35003 


-6.11952 


16.900 


14.259 


12.708 


12.051 


11.804 


M3.2±l.l'= 


N 


-11.3 




Nov 20, 2006 


4 


83.37934 


-6.39085 


18.752 


15.434 


13.776 


13.197 


12.927 


M5.0±0.7'* 


N 


-12.7 




Oct 11, 2011 


5 


83.42722 


-6.05974 


19.915 


16.364 


14.475 


13.786 


13.501 


M4.5±1.0" 


N 


-6.40 




Oct 22, 2011 


6 


83.54100 


-6.38943 


17.206 


14.619 


13.204 


12.576 


12.333 


M3.1±0.5° 


N 


-6.90 




Oct 22, 2011 


7 


83.55022 


-6.28399 


18.653 


15.151 


13.345 


12.711 


12.427 


M5.0±0.7'' 


N 


-9.70 




Oct 11, 2011 


8 


83.56502 


-6.40473 


19.209 


15.796 


13.950 


13.362 


13.053 


M4.0±1.0° 


N 


-5.70 




Oct 22, 2011 


9 


83.63090 


-6.03181 


19.917 


16.230 


14.216 


13.575 


13.287 


M5.1±0.8'= 


N 


-8.70 




Nov 20, 2006 


10 


83.64370 


-6.07540 


19.709 


16.125 


14.116 


13.593 


13.246 


M5.0±2.0° 


N 


-8.70 




Oct 22, 2011 


" J,H,Ks photometry is from 2MASS 
^ Criteria for IR excess are defined in 


Gutermuth et al. 


(2009) and Megeath et al. (2012) 













^ Data from Hectospec spectra. 

Data from IMACS spectra. 
^ Data from [Fang et all j200^') . No spectral-typing error estimate. 
^ Local date at the beginning of the night. 



Table 4 

YSOs with strong emission lines. This table is available only on-line as a machine-readable table. 



ID 


RA 


Dec 


[O I] 


[O I] 


[O I] 


[N II] 


[N II] 


He I 


[S II] 


[S II] 


[O I] 


[O I] 


[Ca II] 


[Ca II] 


[Ca II] 




(J2000) 


(J2000) 


5577A 


6300A 


6363A 


6548A 


6583A 


6678A 


6716A 


673lA 


7773A 


8446A 


8498A 


8542A 


8662A 


6 


83.30896 


-6.19689 






















-0.56 


-0.27 


-0.12 


9 


83.33256 


-6.07309 




-2.0 


-0.37 






















13 


83.35634 


-6.10948 




-5.5 


-1.4 




-0.40 


-1.4 






-0.70 


-1.9 


-6.8 


-6.5 


-5.8 


20 


83.38751 


-6.31044 


-4.7 


-81. 


-27. 




-7.7 


-1.1 


-5.6 


-8.7 


-1.0 


-5.6 


-15. 


-16. 


-13. 


29 


83.42479 


-6.26354 




-4.8 








-2.2 








-6.4 


-5.7 


-6.3 


-5.2 


34 


83.43911 


-6.07385 




-0.49 








-1.8 






-1.0 


-1.8 


-20. 


-22. 


-18. 


38 


83.45954 


-6.36361 




-2.7 


-0.70 




-1.5 




-0.16 


-0.47 












64 


83.57055 


-6.54707 




-4.7 


-1.2 






















89 


83.65066 


-6.09297 




-0.75 


-0.13 






-1.4 






-0.59 


-1.4 


-2.0 


-2.1 


-2.0 


91 


83.65124 


-6.27344 




-1.3 








-2.2 
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Figure 13. Ha equivalent widths (in log scale) of all confirmed 
members vs. spectral type. The red open diamonds are IR-excess 
members and black filled triangles are non IR-excess members. 
The shaded region showrs where the WTTS objects reside using 
CTTS/WTTS classification from White & Basri (2003). Note that 
this classification only applies to stars later than KO. 

widths of the same stars, where the stars with partic- 
ularly low Ha are shown in asterisks. All the low-Ha 
stars shown in the top panel show Li absorption above 
the Pleiades level (|Briceno et al.|[l997l ) and are listed as 
confirmed members. Although it is possible that some 
of these objects are from previous star-forming events in 
the Orion complex and still show Ha and Li, the number 
of low Ha objects is very small and has little effect on 
our overall sample. 

3.4. Extinction 

We use the color excess E(V-I) to estimate the extinc- 
tion toward a star. For sta rs earlier than M4, we use th e 
main-sequence colors from iKenvon fc HartmannI ()1995[ ) ; 
For stars M4 and later, we use the intrinsic colors of 
young disk population described by JjCggctt (1992). We 
then calculate the extinction assuming a standa rd ex- 
tinction law with Ry of 3.1 ()Cardelh et al.l 119890 . The 
E(V-I) metho d has certain l i mitat ions and uncertainties. 
For example, iDa Rio et al.l ()201ClD have shown that ac- 
cretion onto a star makes it look bluer and therefore the 
extinction would be underestimated. 

Figure [T5l shows the distribution of extinction vs. spec- 
tral type . 80% of our spectral-typed members have 
Av< 2. On average, the IR-excess objects (open red 
diamonds) appear to have higher extinction than the 
non-excess objects, which is partially a selection effect 
because we are unable to identify non-excess members 
through Li absorption in high Ay regions where the S/N 
is much lower for the same spectral type. Negative values 
of Ays are non-physical and are due to errors in spectral- 
types, photometry and also uncertainties in their intrin- 
sic colors. Since the spectral- typing errors and photo- 



Figure 14. Top:Ha equivalent widths (in linear scale) of objects 
with Ho equivalent widths less than lA in emission. The dashed 
line is the appr oximate relati on for chromospheric activity in the 
Pleiades (Stauffer et al."1997'). Note that here emission i n H a is 
defined as positive equivalent width, consistent with Figure FlSl but 
different from Table |2] Bottom: Equivalent widths of Li absorp- 
tion for the same objects. Absorption in Li is defined as positive 
equivalent width. The dashed line shows the approximate relation 
for th e maximum L i absorption in the Pleiades for a given spectral 
type IIBriceno et al.i,1997.,') . 

metric errors are larger for the faintest objects, it makes 
sense to see more negative Ays in the latest spectral 
types. 

Figure shows the distribution of extinction of 
spectral-typed members earlier than M5 (solid line) and 
the subset of members earlier than MO (dashed line). 
Stars later than M5 have large uncertainties in their spec- 
tral types and photometry and are excluded. The two 
distributions show that in extincted regions, our sample 
is biased toward more luminous objects because we did 
not apply an extinction cut in our sample selection. 

We now turn our attention to the IR-excess population 
and assess what percentage of this population is highly 
extincted and not in our optical sample. 

We use the 2MASS colors of the IR excess disk objects 
to estimate the extinction and assess the completeness 
of our spectroscopic sample. The top panel of Figure [TTl 
shows J vs. J-H color magnitude diagram of all IR-excess 
objects. The objects that we spectral-typed are shown in 
solid circles (blue) and objects without spectral types are 
shown in open circles (red). Overp lotted are the 3 Myr 
isochrone from iBaraffe et al.l ()1998l) and extinction vec- 
tor from IM© (Ay = 20) and O.IM© (Ay = 10). We use 
the isochrone from Baraffe et al. (1998) since they are 
appropriate for the mass range. We assume an age of 3 
Myr because our optical photometry (Figure [T^ shows 
that LI 641 is approxima tely 3Mrys old (estimated us- 
ing the iSiess et al] j2OOO0 isochrones). and the 3 Myr 
IBaraffe et al.l (|1998D isochrone works the best to match 
up spectral type and J magnitudes. The bottom panel 
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Figure 15. Extinction vs. spectral types for IR excess objects 
(red open diamonds) and non-IR excess objects (black solid tri- 
angles) estimated from E(V-I), assuming the standard extinction 
law from Cardelli et al. (1989) for Ky of 3-1 intrinsic colors 
from|Kenvon & Hartmaim (1995) (for spectra l type s earlier than 
M4) and the young disk population in lLeggetti I I1992I) (for spectral 
types M4 and later). 
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Figure 17. Top: J vs. J-H color magnitude diagram of all IR- 
excess objects. The objects that we spectral-typed are shown in 
solid circles (blue) and objects without spectral types are shown 
in open circles (red ). Overplotted are the 3 Myr isochrone from 
IBarafte et al.l II1998I 1 and extinction vector from IMq (Av = 20) 
and O.IMq {Ay = 10). In the shaded area, our optical spec- 
troscopy sample is about 80% complete. Bottom: J-H vs. H-K 
color-color diagram. The two solid lines are the intrinsic colors 
of main sequence and giant branch ( Bcss ell fc BrettI [T988h . The 
dashed lines and the asterisks show the extinction in of 5 in- 
crements. In both p lots, we use the standard extinction law from 
ICardelli et all iwm for R^ of 3.1. 



Figure 16. Distribution of extinction toward members with spec- 
tral types earlier than M5 (solid line) and members with spectral 
types earlier than MO (dashed line). This illustrates that our sam- 
ple does not have a uniform cutoff in extinction. The extinction 
is estimated from E(V-I) assuming the standard extinction law 
from Cardelli et al. (1989) for Ry of 3.1 and intrinsic colors from 
IKenvon fc Hartmanni (il995f) . Stars later than M5 are not con- 
sidered because of the large uncertainties in spectral-typing and 
photometry. 



shows the J-H vs. H-K color-color diagram. The two 
solid lines are the intrinsic colors of main sequence and 
giant branch (Besse ll fc Br ett 1988,). The dashed lines 
and the asterisks show the extinction in Ay of 5 incre- 
ment s. In both plot s , we u se the standard extinction law 
from lCardeUi et al.l (Il989l) for Ry of 3.1 

From the distribution of the spectral-typed objects, we 
found that our ability to determine spectral types is lim- 
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Figure 18. Same as Figure[T] but the IR-excess stas with spectral 
types are plotted in blue solid symbols and IR excess stars without 
spectral types are plotted in red open symbols. Most of the missing 
spectral types are due to high extinction. 

ited by both the mass (corresponding to the J magni- 
tude) and extinction. For example we can identify and 
spectral type IMq stars in dense regions up to Ay — 8, 
but we can only identify and spectral type O.IMq stars 
up to Ay = 2. There are 234 IR-excess sources in the 
shaded region of the top panel of Figure [17] and more 
than 80% of them (192) are spectral- typed. If we limit 
ourselves to objects with Ay < 2 for 1 to O.IM©, then 
our spectroscopic sample is about 90% (90/110) com- 
plete. We expect similar completeness for non-excess ob- 
jects. Some of the intermediate-mass and massive stars 
are missing as a result of brightness limit of our spec- 
troscopic survey, which we will address in a forthcoming 
paper. Of all the IR-excess objects that lie between the 
O.IMq and IMq line, only 60% of them have Ay > 2, 
which generally agrees with our previous finding that 
only half of the IR excess stars are visible in the opti- 
cal with V< 21. 

These results show that our spectral type distribu- 
tion (Figure 10) is substantially incomplete at the latest 
types. This poses a problem for addressing the shape 
of the IMF at low masses, which we will also address in 
the forthcoming paper. However, from the point of view 
of determining whether the high-mass end of the IMF is 
deficient relative to the low-mass end, missing low-mass 
stars from our sample only means that our results will 
be underestimates of the significance of any depletion of 
high-mass stars. 

Figure [18] shows the distribution of IR-excess objects 
and whether we have their spectral type data. No spec- 
tra were obtained for objects above -6° since they are not 
considered part of L1641. South of -6°, we have spectral 
types for 50% of the disk objects. Most of the missing 
spectral types are due to high extinction. In particular, 
the aggregates of stars between -7.6 - -6 . 8° in declination 
are highly-extincted (jGutermuth et al.ll2011[ ) and there- 
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Figure 19. Top: Distribution of spectroscopically confirmed, 
spectral-typed YSOs (solid line) and all stars with disks (dashed 
line). The star symbols show the number of B0-B4 stars in each 
bin. The low number of spectral-typed stars around Dec = -7.5° - - 
6.8° is a result of high extinction in this declination range. Bottom: 
Cumulative distribution of stars along increasing Declination. The 
horizontal dash-d ot lin e represents a more strict cutoff of L1641 at 
Dec = -6.5° (see [HjJ. 

fore were too faint for spectral typing. 

Figure [TO] shows the distribution in Declination (top) 
and the cumulative distribution (bottom) of the opti- 
cal spectroscopically confirmed members (solid lines) and 
the distribution of all the IR-excess disk stars (dashed 
lines). In low extinction regions, we are able to classify 
most of the class II and class Hi's, and the number of 
spectroscopic members is larger than the IR excess sam- 
ple; whereas in high extinction regions, there are very 
few optical members. For example, the number of opti- 
cal members is very small around Dec = -7.5° - -6.8° due 
to extinction. 

4. DISCUSSION 
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Table 5 

Sample Definition 



Definition 


IR excess disks 


no IR excess 


Total Number 


Spectroscopieally confirmed members 


406"^ 


458 


864 


IR excess sample from Megeath et al. (2011) 


723 




723 


IR excess sample with spectral type 


406= 




406 


IR excess sample without spectral type 


317 




317 


Class Ill's with moderate extinction ^ 




^458 


^458 


All known members ^ 


723 


458 


1181 


Total Number of Class III's'^ 




~900 


~900 


Number ot moderately extincted stars 


406 


;$458 


~ 860 


Total number of stars 


723 


~900 


~1600 



^ Est imat ed number from the V vs. V-I CMD, minus the extrapolated foreground contamination. 
See S im 

^ All known members include all spectroscopieally confirmed members and the IR excess sample 
without spectral type. 
Assuming the class Il/class III ratio is the same for the extincted population. 
This number include 8 objects that are classified as protostars, most of which are in fact 
flat-spectrum objects. 

4.1. Low-Mass Population 

In § 13.21 we discussed our sample of optically con- 
firmed members, including most of the moderately- 
extincted class II's and some class Ill's that were ob- 
served serendipitously. In § 13.41 we discussed the sample 
of highly extincted class II's that are without spectral 
types. In this section, we will estimate the number of 
class III objects expected in L1641 from the CMD and 
compare this number to the number of class Ill's that we 
have already observed. Table [5] gives the definitions of 
the various populations and their numbers. (Note that 
numbers with ^ signs are estimates). 

We noted that most of the IR-excess stars fall in a small 
region on the optical CMD (Figure |4]) between the top 
two parallel lines (described in § 12. ip . We expect that 
the Class III objects have similar ages and distances and 
therefore fall in the same region of the CMD. However, 
not all the non-IR excess objects in this region are mem- 
bers. To estimate the fraction of the non- members in the 
non-IR excess sample, we use the strip one magnitude be- 
low the YSO region as a reference to estimate the number 
of foreground contaminants in the YSO region. Assum- 
ing no foreground extinction and uniform stellar density 
of foreground stars, the distribution of foreground stars 
in a give n magnitude b and is proportional to 10'^'^ (see 
§ 3.6 in iBinnev et al.l [2000'). Therefore, the number of 
foreground stars in the YSO region is predicted to be 
10~ + 10^^'^ times the number of foreground stars in 
the reference strip. 

We restrict our extrapolation to V-I range of 2 - 3.5 
as this is the only portion where we have data for both 
the reference and the YSO region. We identified 514 
stars in the reference strip between V-I of 2 - 3.5 (only 
a very small fraction of which are members). The es- 
timated contamination in the YSO region is therefore 
514 X (lO-^-S + 10-1-2) « 162 stars. There are 496 non- 
IR excess objects in the YSO region between V-I of 2 to 
3.5. We therefore estimate that about 33% of the non- 
excess stars in the YSO strip are non-members. 

The entire YSO strip from V-I of 0.5 to 4 has 772 non- 
excess stars. If we take the same contamination rate of 
33%, we estimate that 772 x (1 - 33%) - 520 of these 
non-excess stars are members. This agrees roughly with 



the number of class III objects (458) that we confirmed 
spectroscopieally and shows that we have identified the 
vast majority of the class III objects subject to our ex- 
tinction limits. 

Approximately half (347 out of 723) of the IR-excess 
objects are identified in the optical photometry. If we 
further assume that the ratio of class II to class Ill's is 
constant even for the more extincted population, then we 
can estimate that there are approximately 458 x 2 ~ 900 
class III objects. The total population of Class II -I- Class 
III stars in LI 641 would then be roughly ~ 1600. 

The above calculation, of course, is a very rough es- 
timate, and depends upon our assumptions. First, we 
assumed that the the foreground stars are all main se- 
quence stars and are uniformly distributed. We have 
also implicitly assumed that the class II to class III ra- 
tio is constant for all spectral types, while the observed 
disk frequency i s a function of spectral type (see e.g. 
iLada et aril2006[ ): however, this is probably a reasonable 
assumption for the majority of our members, as the disk 
frequency seems to decrease most strongly for stars ear- 
lier than M. The assumption that the disk frequency is 
the same for the extincted population can also be prob- 
lematic. The Spitzer/IRAC survey has found that Class 
II objects in L1641 are more distribut ed spatially com- 
pared to the class I's () Allen et al.ll2007l ): similarly, it may 
be that the Class III objects are more distributed and 
have lower extinction systematically than the Class II 
stars. iMegeath et al.l (|2012. in preparation! ) estimated a 
disk fraction of 70 - 80 % in denser regions of L1641 by 
using off fields to subtract out the number of contami- 
nating stars, higher than our the disk fraction we found 
for moderately extincted stars. Therefore, assuming the 
same disk fraction could lead to an overestimate of the 
number of class III stars. 

We can compare the number of low-mass stars in L1641 
to the results from surveys of the ONC. The ONC region 
has been surveyed multiple times both in the IR and in 
the optical. The number of sources found highly depends 
on the area and depth of the specific surveys. Most of the 
IR surveys go much deeper than 2MASS, which makes 
comparing the number of sources m ore difficult. The op- 
tical survey of IDa Rio et al.l (|2010( ) is probably the most 
comparable to ours since they have a 50% V band com- 
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RA (2000) 



Figure 20. Positi ons of known B O - B4 stars near L1641 overlaid 
on the ^3CO map I IBallv et al.lll98 7). The blue solid line outlines 
the Spitzer/IRAC fields. The spectr al type infor mation is taken 
from the spectral type compilation bv \Skiff f201(\) . There are no 
O stars in this region. 

pleteness limit of 20.8. In an area of 34' by 34' they iden- 
tified ~ 1500 stars t hat a re both present in the V and I 
bands. iHillenbrandl (|1997t ) also has a similar survey area 
and found similar number of objects. Therefore, we con- 
clude that L1641 contains comparable number of stars 
as the half degree field centered on the ONC to within 
a factor of 2. This suggests that a direct comparison of 
the stellar IMFs in the two regions can be statistically 
meaningful. 

4.2. High-Mass Population 

With an idea of the number and spatial distribution of 
low-mass stars in L1641, we can now look into the high- 
mass end of the IMF. We search for known early type 
st ars in L1641 in the catalog of spectral types compiled 
bv lSkia Emxh . There are no O stars in the area shown 
and the known early B type (BO - B4) stars are shown i n 
Figure EOl overlaid on the "CO map ()Ballv et al.lll987[ ). 
Table [6] lists the positions, spectral types and photome- 
try of these stars. In the case of spectral classifications 
from multiple sources, the most recent classification is 
used. There are a total of seven early B stars in this 
region. Four of the early B stars (1 - 4 in Table |6]) are 
on the northern tip of the cloud and potentially denote 
the outer regions of the ONC; two other stars (5 and 6) 
are not associated with the cloud if seen in ^'^CO emis- 
sion but could be associated with the lower density gas 
probed by the near-IR (2MASS) extinction map. Since 
we do not have information on the low-mass population 
in their surroundings and they are not clearly associ- 
ated with L1641, we will not consider these two stars 
in the IMF analysis. The other early B star (7) i s th e 
B4 star near (RA,Dec) = (5:42:21.3, -08:08:00) (jRacind 




B-V 

Figure 21. V. vs. B-V CMD of B type stars n ear L1641 (shown 
in Fig ure I20I I . The photometry is taken from I Warren fc Hesserl 
1119771 ). The curved solid line is the ZAMS from Schaller et al. 
1992 (Geneva tracks). Note that B and V photometry are not 
available for all B type stars. The solid line extending from each 
star shows its de-reddened position on the CMD, where the extinc- 
tion is calculated from the E(B-V) color excess. 

[1961 . Figure HH shows the V vs. B-V CMD of these 
early B stars with photometry f rom Warren an d Hesser 
(1977) with the ZAMS from Sch aller et al1[T99l (Geneva 
tracks) overplotted. The dereddened position of each star 
is also shown, where the extinction is calculated from the 
E(B-V) color excess. All the known early B stars are con- 
sistent with the ZAMS at the distance of 420 pc, which 
means they are likely associated with Orion A. 

From Figure [211 we found that all of the known early 
B stars have small extin ctions of Ay < 2 , which is a 
result of the depth of the iWarren fc Hessed (|1977[ ) pho- 
tometry. It is therefore possible that there are embed- 
ded O or early B stars that have not bee n found. How- 
ever, based on the ex isting Spitzer survey ijMegeath et al.l 
12012. in preparation! ) , it is unlikely that we miss an early 
B star. Evolved O or early-B stars would have blown 
away materials nearby in a manner similar to the B4 
star near (RA, Dec) = (5:42:21.3, -08:08:00) and there- 
fore will not have high extinction. The reflection neb- 
ula near the B4 star is the only one seen in our field, 
which supports that there are no evolved early-B stars. 
(Note that there is another reflection nebula, powered 
by a potentially massive star further south, outside of 
the fields of our optical studies.) On the other hand, 
if there were unevolved, highly-extincted early-B stars, 
they would probably heat up the dust nearby and be 
identified as bright IR-excess sources by the Spitzer sur- 
vey. We first examine possible high-mass disk objects 
in the J vs. J-H CMD (top panel of Figure 17). The 
brightest star in the J band is on the upper left corner 
and has a spectral type of Al. A few other stars ap- 
pear to have similar mass and are more extincted, but 
they all have similar masses according to the diagram. 
The only exception is the FU Ori star V883 Ori (with 
J — 9.25 and J-H = 2.49), which is dominated by ac- 
cretion luminosity. We then examine the protostars that 
have large bolometric luminosities. The most luminous 
protostar in L1641 is located at (RA, Dec) = (5:40:27.4,- 
7:27:30) in one of the high-extinction clumps. It has L 
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Table 6 

List of early B-type Stars in L1641 



# 


RA 


Dec 


Name 


Spectral Type'' 




B-V'' 


In Spitzer Field?'^ 




(J2000) 


(J2000) 






(mag) 


(mag) 


Y or N 


1 


83.7542 


-6.00928 


HD 36959 


61.5^* 


5.68 


-0.21 


Y 


2 


83.7612 


-6.00203 


HD 36960 


B0.5V<= 


4.79 


-0.25 


Y 


3 


83.8156 


-6.03275 


HD 37025 


B2/3V' 






Y 


4 


84.1487 


-6.06475 


HD 37209 


B3IVS 


5.72 


-0.22 


Y 


5 


84.6582 


-6.57397 


HD 37481 


B2Vh 


5.96 


-0.22 


N 


6 


85.3434 


-6.93519 


HD 37889 


B2.5Vsn<= 


7.65 


-0.1 


N 


7 


85.5888 


-8.13339 


HD 38023 


B4V' 


8.88 


0.33 


Y 



Most recent spectral type from[Skifi (2010). 
^ Optical photometry from I Warren fc Hesseri II1977I V 

The position on the sky is within the Spitzer field. We use the Spitzer field as a proxy of the main body of the 
cloud. 

The spectral type information comes from Walborn & Fitzpatric3 l lT99Cf) 

The spectral type information comes from Abt & LcvatQ (1977. 1 

f The spectral type information comes from jDcrviz ( 

^ The spectral type information comes from^bt (2008) 

The spectral type information comes from|Zorec ct al.i |[2009l 'l 
' The spectral type information comes from [Racing tl96^ 

~ 490 Lg ijKrvukova et 311120121 ). and according to the 
Siess (2000) isochrones, has a mass less than 7 Mq. 



The lKroupal ()2001l ) IMF can be written as: 



4.3. Comparing the High-Mass IMF of L164-1 to 
Analytical IMFs 

In § 14.21 we discussed the high mass population in 
L1641 from a literature search. In our following analysis, 
we assume that there are no other early B stars in the re- 
gion. All but one of the early B stars are in the northern 
region contiguous with the ONC region. This raises the 
question of how we should define the region in which we 
construct the IMF. We need to adopt a more conserva- 
tive definition and limit our sample to stars from a lower 
declination. This also means we have to compromise for 
a smaller sample size. Figure [19] shows that a large num- 
ber of the confirmed members lie in the northern end of 
the cloud and that our sample size will be considerably 
smaller if we consider only the members further south. 
We need choose carefully a cutoff for our L1641 sample. 

Here we discuss whether this lack of high mass stars 
is statistically significant with the following two ap- 
proaches: 1. considering the entire LI 641 region south 
of -6° and 2. considering only the stars south of -6.5°. 
For comparison, we also consider the region between - 
6° and -6.5°, which is a relatively higher density region 
in L1641 and contains 4 early B stars. Given the sam- 
ple size n and the mass of the most massive star in the 
region Mmax, we can then calculate how unlikely that 
this sample is drawn from analytic IMF models. Here 
we only co nsider the Chabrier (2005) and lKroupal ()2001[ ) 
IMFs. The lChabrie'rl (|2003) IMF gives results that are in 
between the two IMFs discussed here. 

Fi rst we normalize the IMFs between 0.1 and 120Mq. 
The lChabrie^ (pOOSi ) IMF can be written as: 



_ \ = A exp{- 
I = 0.549ylm 



(logm-logO 25)^ 
2x(0.55)^ 
-1.35 



O.IM© < TO < 
m > IM©, 
(1) 



^(to) 



= 2Am 

= A 



1.3 



O.IM© < TO < O.5M0, 
m > 0.5M©, 



(2) 



where the normalization constant A ^ 0.1431. 
The cumulative distribution function is 



I](M) : 



logJ\/ 

^(logm) d log TO 

log(m=O.lM0) 
M 

^(m) dm 

m=O.lM0 



(3) 



where A ^ 0.959 normalizes the total number to unity. 



The probability that a randomly sampled star from 
this IMF is more massive than M^ax is p = 1 — S(Mi„ax)- 
In a population of n stars, the expected number of stars 
more massive than Mmax is -E'(>Mniax) = np. The 
probability that the most massive star is below M,„ax is 
P= (1 -p)". 

Table[7]lists the sample size n, the most massive star in 
the regionMjiiaxand_2 = 1 - S(Mmax ) , -E(> Minax) = np 
for both lChabrieU (l2005l) and lKroupal (l200l IMFs. Note 
that here we consider two sample sizes: spectral- typed 
members and total known members, which is the sum of 
spectral-typed members and the IR-excess members that 
are too ex tincted for optical spectroscopy. Generally, the 
iChabrieii ([2005) IMF gives a higher number of high-mass 
stars and therefore it is less likely to find no high-mass 
stars. 

If we consider the entire L1641 region south of -6°, 
the most massive star is a B0.5V star (corresponding to 
Mmax~ 16M©). We expect to see 2.4 to 4.4 stars more 
massive than the B0.5V star in such a population. If the 
populations are randomly drawn, the probability of not 
finding any star more massive than 16MQis quite low 
(0.012 to 0.091). Therefore, this population has fewer 
-j^^^igh-mass stars compared to the analytic IMF models, 
li\it the result is not very statistically significant. 

If we consider only the region between -6° and -6.5°, 
the most massive star in this region (B0.5V) is compati- 
ble with random sampling of both the analytical IMFs. 
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By limiting our sample to stars south of -6.5° and on 
the cloud (i.e. within the Spitzer snrvcy area), wc can en- 
sure that our sample is not confused with stars from the 
ONC. There are no stars earlier than B4 (corresponding 
to Mmax^ 7M0) in this sample. We expect to see 4.3 
to 9.6 stars more massive than the B4V star in such a 
population. If the populations are randomly drawn, the 
probability of not finding any star more massive than 
7M0is very low (1.4 x 10"^ to 6.4 x 10~^). If we con- 
sider only the 500 members with optical spectral- types, 
the population is incompatible with the analytic IMFs 
to a 2 - 3 (T level; if we further include the 262 IR-excess 
members without optical spectral-type, the population is 
incompatible with the; analytic IMF to a 3 - 4 cr level. 

Even though we currently cannot constrain the total 
number of highly extincted stars, the total population 
will only be larger if we consider highly-extincted class III 
stars. Assuming a complete sample of early B stars and 
the lower limit of low-mass stars, the probabilities we cal- 
culated above are upper limits. Therefore the existence 
of highly-extincted low mass stars will only strengthen 
and not weaken our conclusions. On the other hand, our 
conclusion relies very strongly on the assumption that 
there arc no highly-extinced high mass members, espe- 
cially stars earlier than B4 in the region south of Dec = 
-6.5°. We also want to point out that the expected num- 
ber of high-mass stars varies greatly for different IMF 
models due to their discrepancies in the low-mass end. 
The low-mass IMF is less well-studied and the discrep- 
ancies in the models can be viewed as an indication of 
the uncertainties. 

Since we do not yet have a complete analysis of the 
high-mass and intermediate-mass stars in L1641, we 
choose not to present an HR diagram of the low-mass 
sample alone. We will present the high-to-intermediate 
mass sample in a forthcoming paper, in which we will be 
able to construct a full IMF (down to M3) of LI 641 and 
compare it to that of the ONC. 

5. CONCLUSIONS 

We conducted an optical photometric and spectro- 
scopic survey of LI 641 to test whether IMF varies 
with the environmental densities. Our spectroscopic 
sample consists of IR-excess objects selected from the 
Spitzer /IR AC survey and optical photometry and non- 
excess objects selected from optical photometry. We 
have spectral- typed 406 IR-excess members. We have 
also identified and spectral-typed 458 non-excess YSOs 
and 98 probable members through Li absorption and Ha 
emission. Our sample is limited by extinction and stellar 
mass in a complicated way, but overall we are able to 
spectral-type about 90% of class II and class III YSOs 
between 0.1 and 1 Mq with Ay ^ 2. The number of 
spectroscopically confirmed class III objects is also con- 
sistent with the number of class III objects estimated 
from extrapolation of the V vs. V-I color magnitude di- 
agram. 

Assuming the class Ill/class II ratio is the same for the 
more extincted population, we estimate the total number 
of class II and class Ill's in L1641 to be around 1600, 
including 723 class II objects from the Spitzer survey 
and ~900 class Ill's. 

The total number of stars in L1641 is comparable 
within a factor of two to the number of stars used to 



construct the IMF in the ONC, even though the number 
depends on how our L1641 sample is defined. The op- 
tical photometry shows that the ONC and LI 641 are of 
similar age, with L1641 being slightly older by ~1 Myrs, 
indicating that any late-0 to early-B star should still be 
on the main sequence. Assuming that we know all the 
early B stars in L1641, we compare the high-mass IMF to 
the standard analytical models. Compared to the stan- 
dard models of the IMF, L1641 is deficient in O and early 
B stars to a 3-4cr significance level. We will make further 
improvements to our sample to better constrain the ex- 
tincted population and compile a complete sample of B 
stars in a forthcoming paper, and the we will be able to 
make a direct comparison with the ONC. 

We thank the referee for helpful suggestions to improve 
the manuscript. W.H. and L.H. acknowledge the sup- 
port of NSF grant AST-0807305 and the Origins grant 
NNX08AI39G. This paper uses data obtained at the 
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stitution and the University of Arizona and the 6.5 meter 
Magellan Telescopes located at Las Campanas Observa- 
tory, Chile. This paper uses data products produced by 
the OIR Telescope Data Center, supported by the Smith- 
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Table 7 

Samples by Region and IMF Probabilities 



Kegion 


n (Sample Size) 




p=l- S(Mmax)'' 


E{> Mmax)"^ 


P(no stars greater than Mmax)'^ 


L1641 


S])cctral-typc<i members 864 


B0.5V 


ChabrierOS 3.74 X 10-^ 


3.2 


0.039 


Dec < -6° 




(KiAI. ) 


Kr(.m|)a01 2.77 X 10 


2.1 


0.0!)1 




Total Known Members^llSl 


B0.5V 


Cliabricr05 3.74 X lO"'^ 


4.4 


0.012 






(I6M0) 


KroupaOl 2.77 x 10"^ 


3.3 


0.038 




Spectral-typed members 364 


B0.5V 


Chabrier05 3.74 x lO"'' 


1.4 


0.26 


Dec = -6° to -6.5° 




(16Mq) 


KroupaOl 2.77 X IQ-^ 


1.0 


0.36 




Total Known Members" 419 


B0.5V 


ChabrierOS 3.74 x lO"'^ 


1.6 


0.20 






(I6M0) 


KroupaOl 2.77 x IQ-^ 


1.2 


0.31 




Spectral-typed members 500 


BIV 


ChabrierOS 1.26 x 10"^ 


6.3 


1.8 X 10"-^ 


Dec < -6.5° 




(7M0) 


KroupaOl 8.55 x IQ-^ 


4.3 


1.4 X 10-2 




Total Known Members" 762 


B4V 


Chabrier05 1.26 x 10"^ 


9.6 


6.4 X 10-° 






(7M0) 


KroupaOl 8.55 x IQ-^ 


6.5 


1.4 X 10-3 



Most massive star observed in this region. 
^ The probability that a random star drawn from the IMF is more massive than Mmax ■ 
Expected number of stars more massive than Mmax- E{> Mmax) = np 

Drawing n stars from the given IMF, the probability of getting no stars more massive than Mmax- P = (1-p)'^ 
" Total known members include spectral-typed members and IR excess members without spectral type. 
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